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ISM ENRICHMENT AND LOCAL POLLUTION IN DWARF 

GALAXIES 

D. Kunth 1 and V. Lebouteiller 2 



Abstract. The fate of metals after they are released in starburst episodes 
is still unclear. What phases of the interstellar medium are involved, 
in which timescales? Evidence has grown over the past few years that 
the neutral phase of blue compact dwarf (BCD) galaxies may be metal- 
deficient as compared to the ionized gas of their H II regions. These 
results have strong implications for our understanding of the chemical 
evolution of galaxies. We review here the main results and the main 
caveats in the abundance determination from far-UV absorption-lines. 
We also discuss possible scenarios concerning the journey of metals into 
the interstellar medium, or even their ejection from the galaxy into the 
intergalactic medium. 



1 Introduction 



p<| ■ Within the chemical downsizing scenario, in which massive galaxies are the first 



o 



to form stars at a high rate (Cen & Ostriker 1999), dwarf galaxies can remain 
chemically unevolved since their formation. It is thus possible that some dwarf 
galaxies in the nearby Universe still contain pristine gas (Kunth & Sargent 1986; 
Kunth et al. 1994). 

The family of Blue Compact Dwarf (BCD) galaxies represents a major testbed 

■ as it includes the most-metal poor star-forming galaxies known (Kunth & Ostlin 
2000). Most of these galaxies contain young (3-10 Myr) super stellar clusters that 

■ dominate the interstellar medium (ISM) enrichment in a-elements. Old stellar 
populations have been detected in BCDs (e.g., Doublier et al. 2001; Aloisi et al. 
1999; Aloisi et al. 2001; Ostlin 2000; Vanzi et al. 2000; Hopp et al. 2001), sug- 
gesting that although "chemically young", these objects are not genuinely young. 
Nevertheless, only 22% of the total stellar mass in I Zw 18 could be contributed by 
older stars which in itself makes this object peculiar per se (Hunt et al. 2003). 
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Although evidence of local enrichment has already been observed (Walsh & 
Roy 1993), it is usually associated with the presence of Wolf-Rayet (WR) stars 
(Walsh & Roy 1989; Thuan et al. 1996; Lopez-Sanchez et al. 2007). Brinchmann et 
al. (2008) confirmed this assumption by analyzing 570 WR galaxies showing larger 
N/O abundance ratios as compared to normal galaxies. In order to probe the ISM 
enrichment by starburst episodes, dispersion and mixing timescales have to be 
better understood. For example, a typical lOMyr burst involving 2.5 x 10 6 M Q 
of stars at a metallicity of 0.2 Zq yields about 10 4 Mq into the ISM (using yields 
of Meynet & Maeder 2002). Considering a mass of ionized gas of 2.5 x 1O 5 M0, 
an instantaneous local enrichment of the newly produced metals would lead to 
an unrealistic metallicity of 10 Z Q ! It is clear that observing the ionized gas only 
(through the optical-infrared emission-lines) does not provide the complete picture 
of the journey of heavy elements into the ISM. The neutral envelope notably plays 
a fundamental role. 

The role of the neutral gas upon the chemical evolution of galaxies is twofold. 
First, it acts a massive reservoir for the molecular hydrogen which fuels star- 
formation. Accretion episodes provide theoretically enough gas to explain the 
star- formation rate observed in certain galaxies (e.g., Fraternali & Binney 2008). 
The importance of accretion is likely enhanced in dwarf galaxies characterized by 
a lack of rotation and inefficient cooling flows. Furthermore, as star-formation 
occurs within giant H n regions, newly produced metals can be driven to high 
galactic latitudes into the H I halo through the so-called galactic fountains. The 
dispersion spatial and time scales as well as the mixing time scale are however 
relatively unknown in dwarf galaxies. 

The sample of BCDs for which the comparison was made between chemical 
abundances in the H I envelope and in the ionized gas of their H II regions is 
presented in TableQ] It includes the 2 lowest-metallicity star-forming galaxies 
known I Zw 18 (Aloisi et al. 2003; Lecavelier et al. 2004) and SBS 0335-052 (Thuan 
et al. 2005), as well as more metal-rich objects, NGC 1705 (Heckman et al. 2001), 
Mrk59 (Thuan et al. 2002), IZw36 (Lebouteiller et al. 2004), NGC 625 (Cannon 
et al. 2004), and Pox 36 (Lebouteiller et al. 2009). In addition, the giant H n 
region NGC 604 in M33 was investigated by Lebouteiller et al. (2006). 

We first explain the strategy for determining abundances in the neutral gas of 
galaxies (Sect.[2|). One major challenges concerns the absorption-line contamina- 
tion by stellar photospheres for the bursts older than lOMyr (Sect.[3|). We then 
review the results in Sect. [Hand investigate the effect of depletion on dust grains 
in Sect.O Finally, we gather some considerations on the enrichment of the neutral 
gas in Sect. El 

2 Determining abundances in the neutral gas 

Metals from the neutral gas can be observed through resonant lines in the far- 
ultraviolet (FUV) . The most massive stars present in BCDs provide a strong FUV 
continuum, and are used as lighthouses to probe the foreground interstellar gas. 
The UV photons cross respectively the stellar atmospheres (absorption of H I and 
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Table 1. Sample of star-forming objects whose neutral gas was investigated with FUSE. 
GHIIR stands for giant H n region. [O/H] is defined as log(0/H) — log(O/H) . We use 
the solar reference of Asplund et al. (2005). 



Object 


Type 


[O/H] (ionized gas) 


IZwl8 


BCD 


-1.45 


SBS 0335-052 


BCD 


-1.36 


I Zw36 


BCD 


-0.89 


Mark 59 


BCD 


-0.67 


Pox 36 


BCD 


-0.61 


NGC625 


BCD 


-0.47 


NGC 1705 


BCD 


-0.45 


NGC 604 


GHIIR 


0.00 



highly ionized species such as C in, C iv, N v, ...), the photoionized gas (C II, 
C in, S in, ...), the neutral ISM (H I, N I, O I, Si II, P II, Ar I, Fe II, ...), and 
the Milky Way ISM. Provided the radial velocity of the object is large enough to 
distinguish it from the Milky Way component, the neutral gas is easily probed in 
BCDs, which display large amounts of H I gas (Thuan & Martin 1981). 

The FUSE telescope (Moos et al. 2000) gave access to the range 900-1200A at 
a spectral resolution of 15 000, allowing for the first time precise determination of 
metal abundances in the neutral gas, following the pioneering work of Kunth et al. 
(1994) with the HST/GHRS. Although several neutral gas tracers are available, 
their use to derive chemical abundances can be challenging for many reasons: 

• The H I lines show prominent damping wings from which the H I column 
density is easily determined. All the lines from the Lyman series fall within 
the FUSE range except Lya. The intrinsically large line width makes it often 
difficult to separate from the Milky Way H I lines. Moreover, stellar atmo- 
spheres can contaminate significantly the absorption profile (more details in 
Sect.©. 

• O I is an important tracer of the a-elements produced in massive stars. The 
ionization potential of O I is very close to that of H I (Table[2|) and the charge 
exchange between O II and H I is very efficient so that the coupling of the 
oxygen and hydrogen ionization fractions is quite strong. Yet, practically, 
O I lines are often saturated due to strong oscillator strengths, making the 
column density measurement uncertain. 

• N is produced by intermediate-mass to massive stars. The presence of N I 
in the neutral regions is well established (Tabled, it can be shown that the 
applied correction to account for the possible presence of ionized N in the 
neutral gas is only 0.05 dex (Sofia & Jenkins 1998). Finally, little depletion 
on dust is expected. An example of N I line detection is shown in Fig. [I] 
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Table 2. Ionization potentials (H I = 13.60 eV). The species observed with FUSE are in 
bold. 



I 


ii 




in 




O I = 13.62 eV 


ii = 


= 35.12eV 






N I = 14.53 eV 


N ii = 


= 29.60 eV 






Si i = 8.15 eV 


Si ii 


— 16.34 eV 


Si in = 


34.49 eV 


P i = 10.49 eV 


P II : 


= 19.72 eV 


P in 


30.18 eV 


Ar I = 15.76 eV 


Ar II 


= 27.63 eV 






Fe i = 7.87 eV 


Fe II 


— 16.18eV 


It hi = 


= 30.65 eV 



• Si II is the main ionization stage of Si in the neutral gas due to the very low 
ionization potential of Si I (TableEJ . Still, ionization correction are necessary 
due to the presence of Si II in the ionized gas. Silicon is depleted on grains, 
though usually less than iron (e.g., Sofia et al. 1994; Savage & Sembach 
1996). 

• P II and O I column densities are seen to follow each other in various en- 
vironments (Lebouteiller et al. 2005; but see Jenkins 2009). However, the 
ionization potential of P n suggests that ionization corrections might be 
needed due to the presence of P II in the ionized gas. 

• Ar is produced in massive stars and should not be depleted in a low-density 
H I cloud (Jenkins et al. 2000). Ionization correction from low-density, 
partly-ionized regions of the ISM can be required because of the relatively 
large photoionization cross section of Ar I (Sofia & Jenkins 1998). 

• Fe II is the main ionization stage in the neutral gas. Fe is significantly 
depleted on dust grains, but a variety of Fe n lines with different oscillator 
strengths is available within the FUSE range, making the observation of 
weak (not saturated) lines possible. 

Chemical abundances in the neutral gas are simply calculated using the col- 
umn density ratio of metals to H I, with some ionization corrections when needed 
(concerning mostly Ar I, Si II, P n). We refer to the papers on the individual 
objects for more details. 

3 H i contamination from the stellar atmosphere 

The FUV spectrum of most of the BCDs and giant H II regions studied so far 
with FUSE is dominated by O stars. This is indicative of a starburst age younger 
than < lOMyr (Robert et al. 2003). Because of the high temperature, hydrogen 
is mostly ionized, and only weak H I photospheric lines can be observed which 
contribute to the — already saturated — core of the interstellar H I absorption 
line. Hence FUV spectra toward young starbursts allow a precise determination 
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Fig. 1. N I triplet at 1134 A toward IZwl8 (figure taken from Lecavelier et at 2004). 

of the interstellar H I column density, iV(H i), by using the damping wings. It 
must be noted however, that extremely young starburst episodes (~ 3 Myr) show 
a prominent O VI P-Cygni profile which complicates the profile fitting of the Ly/3 
line (Lebouteiller et al. 2006). 

Old starbursts (> 10 Myr) are dominated by B stars. Such stars have overall a 
more complex FUV spectrum than O stars because the ionization degree is lower. 
For this reason, photospheric H I lines becomes prominent (Valls-Gabaud 1993). 
As an illustration, the H I lines from the Lyman series in Mrk 59 and Pox 36 show 
"V-shaped" wings typical from B stars photospheres. While Thuan et al. (2002) 
circumvented the problem by considering an artificial continuum to fit the profile 
of the H I lines in Mark 59, it is also possible to model the stellar absorption. 

In their study of Pox 36, Lebouteiller et al. (2009) used the TLUSTY models 
for O and B stars (Lanz & Hubeny 2003; Lanz & Hubeny 2007) to reproduce the 
synthetic spectral continuum. The stellar spectrum was modeled with a single 
stellar population. It appears that the wings of the stellar H I lines are strongly 
dependent on the star temperature. The best constraint is given by the high-order 
line of the Lyman series since the interstellar contribution is saturated with the 
absence of damping wings (Fig. [2]). 

On the other hand, since low-order Lyman lines (Lya, Ly/3) are characterized 
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Fig. 2. H I Lye line toward Pox 36 (from Lebouteiller et al. 2009). The "V-shaped" wings 
are typical of stellar photosphere absorption and their shape is strongly dependent on 
the star temperature. The core of the line is saturated and originates from the narrow 
interstellar H I. In the case of Pox 36, the dominant stellar population has T ~ 30 000 K 
(B0). 

by strong interstellar damping wings, interstellar H I dominates the global shape. 
For instance, Ly/3 can be used to constrain the interstellar H I column density if 
the stellar temperature is > 30 000 K (Fig. [3]). 

As a summary, the profile of Lya is a good constrain on iV(H i) regardless 
of the burst age while the profile of Ly/3 can be used only for a dominant stellar 
population warmer than 30 000 K. Higher-order Lyman line profiles provide only 
little constraints to 7V(H i) if the signal-to-noise ratio is less than about ~ 50. 

4 Results 

Results are shown in Fig. HI where the abundances of N, O, Ar, and Fe in the 
neutral gas of BCDs are compared to the abundances in the ionized gas of the H n 
regions. The same dataset was used by Aloisi et al. (2003) and Lecavelier et al. 
(2004) to investigate IZwl8 (data points (1) and (1*) respectively). Lecavelier et 
al. find O i/H I— — 4.7^g which is consistent with the O/H ratio observed —0.6 
in the H n regions (all uncertainties are 2-a) while Aloisi et al. report a significantly 
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Fig. 3. Complicated figure showing the H I Lyman line profiles (stellar and stel- 
lar+interstellar) as a function of the star temperature (horizontal), of the Lyman or- 
der (vertical), and of the interstellar H I column density (line styles) (figure taken from 
Lebouteiller et al. 2009). For hot stars, the profile of Lya and Ly/3 is dominated by the 
interstellar component. For stars colder than ~ 30 000 K, only Lya can be used for a 
reliable interstellar H I column density determination, as the stellar absorption starts to 
dominate the line profile. 



different value with O i/H i= —5.4 ±0.3, a discrepancy with the former result that 
awaits for clarification. Hereafter, we define the underabundance of an element X 
in the neutral has as: 



5 H i(X) = [X/H] HII -[X/H] 



(4.1) 
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where [X/H]h n is the abundance in the ionized gas of the H II regions and [X/H]h i 
is the abundance in the neutral gas. 

Examination of [Ar/H] and [N/H] in Fig. 0] reveals the possible existence of a 
plateau for the lowest- metallicity BCDs. The data on [O/H] are also consistent 
with the existence of a metallicity plateau located at 12 + log (O/H) ~ 7.0 ([O/H] 
= —1.7), although the error bars are considerably larger. Part of the dispersion 
of [O /H] is probably due to saturation effects which could not be avoided in some 
objects. Some of the dispersion of [N/H] and [Ar/H] is probably due to ionization 
corrections. 

Two fundamental results can be drawn from Fig. 01 (1) the abundances in the 
neutral gas show that this phase has been already enriched in heavy elements, and 
(2) the abundances in the neutral gas are systematically equal or lower than those 
in the ionized gas (<5hi(X) > 0). These results suggest that the neutral envelope 
of BCDs is participating in the dispersion and mixing of heavy elements produced 
in massive stars. However, the neutral phase seems to be less chemically evolved 
than the ionized gas associated with the current star-formation episode. 

The interpretation of abundances in BCDs has been limited so far by the lack 
of knowledge about the origin of the neutral gas whose absorption is detected in 
the FUV. While absorption-lines can arise in the cold and warm neutral medium, 
it must be kept in mind that there is a bias in the FUV to observe dust-free lines 
of sight toward the most massive stars. The direct comparison with the H I gas 
detected via the 21 cm emission in the BCDs is also hampered by the fact that FUV 
observations only probe gas in front of stellar clusters while radio observations (in 
emission) probe the whole system. We discuss in the following the properties of 
the FUV lines of sight by the means of depletion calculation. 

5 Depletion on dust grains 

Depletion concerns mostly Si, Fe, and to a minor extent, P. The depletion strength 
in the neutral envelope of BCDs is largely unknown, mostly because of the lack of 
information on the FUV lines of sight and because of the low metal content of the 
gas. This is an important parameter to account for to estimate the fraction of the 
metal underabundance in the neutral gas genuinely due to the chemical evolution. 
The observed abundance of an element X can be expressed as: 

[X/H] « log(Z) + [X/H] Q + A X F«, (5.1) 

where [X/H] is the abundance respect to the solar value, i.e., log(X/H)— log(X/H)©. 
Z is the intrinsic metallicity of the gas in solar units, [X/H]o corresponds to the 
initial (minimal) depletion, Ax (< 0) is the propensity of a given element to be 
depleted and F* is the sightline collective depletion strength, following the defi- 
nitions of Jenkins (2009). It must be noted that [X/H]o depends on metallicity, 
being likely closer to in metal-poor environments because of the lower dust-to-gas 
ratio (Hirashita et al. 2002; Lisenfeld & Ferrara 1998). 

For simplicity, we make the assumption that the metallicity of the neutral gas 
is solar log(Z) = and try to find the depletion strength F* agreeing with the 
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Fig. 4. Abundances of N, O, and Ar in the neutral gas and in the ionized gas (figure 
taken from Lebouteiller et al. 2009). The dashed line indicates the 1:1 ratio. Labels: 
(1*) and (1) IZwl8 (see text for details), (2) SBS 0335-052, (3) IZw36, (4) Mark 59, (5) 
Pox 36, (6) NGC625, (7) NGC 1705, (8) NGC604/M33. 
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Fig. 5. Depletion strength in the neutral gas assuming that it is the parameter dominat- 
ing the observed abundances in the neutral gas. This is an upper limit since we assumed 
a solar metalhcity and the neutral gas is intrinsically metal-poor (see text). 



observations. This is equivalent to assuming that #hi is dominated by depletion 
on dust grains rather than a different chemical enrichment. Figure[5] shows that 
F 1 * in BCDs would then be comprised between 0.2 and 1.1 (based on the use of 
the most refractory element, Fe). This is a strong upper limit because we assumed 
a solar metalhcity. In practice, to recover the observed metallicity in the ionized 
gas, F, would be essentially null. The fact that F* = would imply that the lines 
of sight are essentially dust free, intersecting UV transparent diffuse clouds. This 
is of course compatible with the selection effect of lines of sight in the sense that 
the most opaque clouds are invisible in the UV. 

If F* is indeed null, underabundances in the neutral gas are not due to depletion 
on dust grains, as it is already suggested in Fig. [4] by the study of N and Ar (only 
little depleted). The observed abundances of N, O, and Ar can only be explained 
by a genuine chemical underabundance in the neutral gas. 

The situation for Fe requires further interpretation. Fe is known to deplete 
significantly in the ionized gas, even in metal-poor objects. In the Milky Way, it 
is thought that about 2/3 of the iron in the ISM was produced by SNe type la 
involving low-mass stars. Rodriguez & Esteban (2004) and Izotov et al. (2006) 
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Fig. 6. Abundances of Fe in the neutral and ionized gas. See Fig[4]the plot description. 
Grey points correspond to the [Fe/H] in the ionized gas (observed values) while black 
dots include a correction to account for depletion in the ionized gas (figure taken from 
Lebouteiller et al. 2009). 



observed that the [O/Fe] ratio increases with metallicity in emission-line galaxies, 
being consistent with depletion of Fe on dust grains in the ionized gas, and with 
iron production in SN type II in those galaxies. Hence we could be in a situation 
where Fe is less depleted in the neutral gas — probed along dust-free lines of sight 
— than in the ionized gas. Correcting for depletion of Fe in the ionized gas leads 
to similar abundance patterns to the other elements, i.e., an underabundance in 
the neutral gas (see Fig. [6] and Lebouteiller et al. 2009). 



6 Discussion 



We now assume that all the problems related to the chemical abundance deter- 
minations in the neutral gas are solved. These include ionization correction, line 
saturation effects, depletion in dust grains, contamination by stellar photospheres 
(for H i). Can the results be understood in a chemical evolution scenario involving 
both the ionized gas and the neutral gas? 
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6.1 Abundance discontinuity 

Considering the mass of metals released by a starburst, several hypotheses come 
into play, as explained in Lebouteiller et al. (2008): 

• All the metals mix in the H n region alone, leading to unrealistically high 
abundances (see introduction). 

• All the metals mix in the H I envelope. Given the large mass of atomic hydro- 
gen, the metallicity of the neutral gas would increase by ~ 5% (calculation 
made for Pox 36). 

• Only a fraction of the metals mixes within the H n gas, and the remaining 
mix within the H I gas. It is possible that the H I close to the clusters is 
enriched first. Lebouteiller et al. (2009) calculated that if metals mix with 
10% of the H I in Pox 36, the enriched H I would acquire the metallicity 
observed in the ionized gas of the H II regions. 

The amount of metals released by a single starburst episode cannot be mixed 
entirely within the H II gas since it would lead to unrealistic abundances. The 
minimal enrichment therefore applies to only a fraction of the expelled metals. A 
possible scenario might thus include a local enrichment in the H n region of a — 
small — fraction of the metals released by the current star-formation episode while 
most of the metals mix in the H I envelope. In order to reproduce the observed 
abundance discontinuity between the 2 phases, the fraction of metals mixing in 
the neutral gas has to satisfy to the condition: 

fn = 1 I 10 fl„, M(HH)- > 
1 + lu M(HI) 

where M (H n) and M(H i) are the masses of the ionized gas and neutral gas respec- 
tively (see Lebouteiller et al. 2009) and where 5hi is the abundance discontinuity 
between the ionized gas and neutral gas as defined in Sect.HJ 

In Pox36, dm m 10 translates into /„ of about 1%. In IZwl8, <5hi ~ 
(following the result of Lecavelier et al. 2004) translates into a fraction f n of 5%. 
The large difference in behavior between /„ and i5hi in these 2 galaxies is due to 
the fact that the ratio of ionized gas mass over neutral gas mass is much larger in 
I Zw 18 than in Pox 36. 

6.2 Metallicity threshold? 

The neutral gas of BCDs has already been enriched with metals up to an amount 
of > 1/50 Zq. Could this value be a metallicity threshold? The hypothesis of 
a minimal metallicity has been already suggested by Kunth & Sargent (1986) to 
explain the difficulty to find extremely-metal poor star-forming objects. In these 
objects, the current star-formation episode could pollute the ISM by heavy ele- 
ments released by stellar winds and supernovas. It is interesting to notice that the 
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apparent threshold of > 1/50 Z in the neutral gas of BCDs is also the metallicity 
of the most metal-poor galaxy known, IZwl8. It is conceivable that the results 
are related. However, they imply considerably different (both spatial and time) 
scales for the dispersion and for the mixing. 

Interestingly, Telfer et al. (2002) found that the intergalactic medium (IGM) in 
the redshift range 1.6 < z < 2.9 has a metallicity 12 + log(0/H) between 6.7 and 
7.6, which includes the floor metallicity observed in BCDs. Hence it appears that 
the neutral and ionized ISM of the lowest metallicity galaxies could well be charac- 
terized by a default enrichment rather than by previous star-formation episodes. 
In that view, BCDs could have started to form stars only "recently" (around 
~ 8 — 9Gyrs ago), from an H I cloud already at the IGM metallicity. This has 
important consequences on the number of starburst episodes that occurred until 
now. For example, assuming /„ = 5% in IZwl8 (to reproduce the result (5hi = 0, 
Sect. 16. 1|) . 3 bursts are required to explain the observed metallicity if starting from 
a null metallicity (calculation made using yields from Meynet & Maeder (2002) 
and assuming interphase mixing between 2 consecutive bursts). However, a single 
burst is required if starting from Z < 1/500. Further observations and analysis 
of BCDs with metallicities lower than ~ 1/8 Zq will be needed to confirm the 
existence of this plateau. 

6.3 Metals from star-forming dwarfs: retention or ejection? 

Since a local enrichment of 100% of the products released by a starburst would 
lead to unrealistically high abundances, 2 possibilities remain: (1) metals disperse 
and mix in another phase (Sect. lfTTj) . or (2) most of the metals are ejected into 
the IGM. We now discuss the parameters playing a role in the ejection/retention 
of metals. 

As they end their lives massive stars explode as a supernovae. The energy 
output from a SN is over a short period, comparable to that of a whole galaxy. In 
a galaxy with a high local star formation rate, the collective action of supernova? 
may lead to a galactic superwind, which may cause loss of gas. Stellar winds can 
also contribute to the energetics of the ISM at the very early stage of a starburst 
(Leitherer et al. 1992). The relative importance of stellar winds compared to SNe 
increases with metallicity. A continuous wind proportional to the star formation 
rate has been applied in models predicting the evolution of starburst galaxies. But 
since different elements are produced on different timescales, it has been proposed 
that only certain elements arc lost (or in different proportions) hence reducing 
the effective net yield of those metals as compared to a simple chemical evolution 
model (Matteucci & Chiosi 1983; Edmunds 1990). The SNe involved in such a 
wind are likely to be of type II because type la SNe explode in isolation and will 
less likely trigger chimneys from which metals can be ejected out of the plane of 
a galaxy. In this framework O and part of Fe are lost while He and N (largely 
produced by intermediate stars) are not. This would result in a cosmic dispersion 
in element ratios such as N/O between galaxies that have experienced mass loss 
and those that have not. In a dwarf galaxy which has a weaker gravitational 
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potential, these effects may result in gas loss from the galaxy unless as we argue 
below the presence of a low H I density halo acts as a barrier. 

Galactic winds have been observationally investigated in dwarf galaxies (e.g., 
Marlowe et al. 1995; Martin(1996; Martin(1998) and more recently with the advent 
of the Chandra satellite. Lcqucux et al. (1995), Kunth et al. (1998), and Mas- 
Hesse et al. (2003) have shown that the escape of the Lya photons in star-forming 
galaxies strongly depends on the dynamical properties of their interstellar medium. 
The Lyman alpha profile in the BCG Haro 2 indicates a superwind of at least 
200 km s -1 , carrying a mass of « 1O 7 M , which can be independently traced 
from the Ha component (Legrand et al. 1997). However, high speed winds do 
not necessary carry a lot of mass. Martin (1996) argues that a bubble seen in 
IZw 18 (see also Petrosian et al. 1997) will ultimately blow-out together with its 
hot gas component. Although little is known about the interactions between the 
evolving supernova remnants, massive stellar bubbles and the ISM it is possible 
that an outflow takes the fresh metals with it and in some cases leaves a galaxy 
totally cleaned of gas. This scenario clearly contradicts the self enrichment picture 
— unless a fraction of hot gas rapidly cools down — . But since we see some hot 
gas outside of the H n regions the question remains whether this gas will leave a 
galaxy or simply stay around in the halo? 

Model calculations developed by Silich & Tenorio-Tagle (1998) predict that su- 
perbubbles in amorphous dwarf galaxies must have already undergone blowout and 
are presently evolving into an extended low-density halo. This should inhibit the 
loss of the swept-up and processed matter into the IGM. Recent Chandra X-rays 
observations of young starbursts indicate some possible metal losses from disks 
(Martin et al. 2002) but not in the case of NGC 4449 with an extended H I halo of 
around 40kpc (Summers et al. 2003). In a starburst galaxy, newly processed ele- 
ments are produced within a region of ~ 100 pc size. During the supernova phase 
the continuous energy input rate from coeval starbursts or continuous star-forming 
episodes, maintains the temperature of the ejected matter above the recombination 
limit (T ~ 10 6 K) allowing superbubbles to reach dimensions in excess of lkpc. 

The mechanical energy released during a starburst episode accelerates the in- 
terstellar medium gas and generates gas flows. The properties and evolution of 
these flows ultimately determine the fate of the newly formed metals and the man- 
ner they mix with the original interstellar medium. The presence of outflows may 
indicate, at first, that supernova products and even the whole of the interstellar 
medium may be easily ejected from the host dwarf systems, causing the contamina- 
tion of the intra-cluster medium (Dckel & Silk 1986; De Young & Heckmanl994). 
This type of assumption is currently blindly used by cosmologists in their model 
calculations. However the indisputable presence of metals in galaxies implies that 
supernova products are not completely lost in all cases (Silich & Tenorio-Tagle 
2001). Legrand et al. (2001) have compared Silich & Tenorio-Tagle (2001) theo- 
retical estimates with some well-studied starburst galaxies. They have worked out 
three different possible star formation history scenarios that assume either a very 
young coeval starburst or extended phases of star formation of 40Myr or 14Gyr 
and inferred the expected energy input rate for each galaxy using the Ha lumi- 
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Fig. 7. Log of the critical mechanical luminosity (left-side axis) and mass of the star clus- 
ter Msb (right-side axis), required to eject matter from galaxies as a function of Mism 
(figure taken from Legrand et al. 2001). Lower limit estimates are shown for galaxies 
with extreme ISM density distributions: flattened disks (lower two lines), and spheri- 
cal galaxies without rotation (upper lines), for two values of the intergalactic pressure 
PiGvi/k = 1 cm -3 K (solid lines) and PiGm/k — 100 cm -3 K (dashed lines). 



nosity and/or the observed metallicity. Values of the derived mechanical energy 
injection rate in the three considered cases were compared with the hydrodynam- 
ical models predictions of Mac Low & Ferrara (1999) and Silich & Tenorio-Taglc 
(2001). Detailed calculations and accompanying figures are given in Legrand et al. 
(2001); we only present in Fig.[7]the resulting plot for the extended phases of star 
formation of 40 Myr that might be the most realistic case for starburst galaxies. 

The net result (see Fig. [7]) is that all galaxies lie above the lower limit first 
derived by Mac Low & Ferrara (1999) for the ejection of metals out of flattened 
disk-like ISM density distributions energized while most are below the limit for 
the low density halo picture. Thus the mass of the extended low density halo 
efficiently acts as the barrier against metal ejection into the IGM. 

Disk-like models clearly require less energy to eject their metals into the inter- 
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galactic medium because the amount of blown out interstellar gas that can open 
a channel into the intergalactic medium is much smaller than in the spherically 
symmetric limit, where all of the metal-enriched ISM has to be accelerated to 
reach the galaxy boundary. Predicted haloes, despite acting as a barrier to the 
loss of the new metals, have rather low densities (< rihaio ~ lCT 3 cm~3) and thus 
have a long recombination time (t rcc = l/(cmh a io); where a is the recombination 
coefficient) that can easily exceed the lifetime of the H II region (foil = 10 7 yr) 
developed by the starburst. In such a case, these haloes may remain undetected at 
radio and optical frequencies until large volumes are collected into the expanding 
supershells. 

6.4 Are BCDs special? 

Bowen et al. (2005) found similar abundances in the neutral and ionized gas of 
the dwarf spiral galaxy SBS 1543+593. This results poses a fundamental ques- 
tion: is the metal deficiency of the H I gas specific to BCDs? Rotation in spirals 
provide an efficient mixing throughout the galaxy. Furthermore, star-formation 
keeps occurring along the spiral arms. In opposition, BCDs have weak velocity 
field, and star- formation usually occurs in one or several knots (e.g., Papaderos 
et al. 1996). The different mode of star- formation could be responsible for sub- 
stantial differences in the dispersal and mixing mechanisms. Lebouteiller et al. 
(2006) found that the neutral phase toward NGC 604 in the spiral galaxy M 33 
is more metal-poor than the ionized gas of the H n region, i.e., being consistent 
with results in BCDs, but the contamination by high-velocity clouds along the 
lines of sight similar to the ones observed around the Milky Way complicates the 
interpretation. It seems necessary to study the metal abundance in the neutral 
gas of dwarf irregular galaxies that are not BCDs. Studies of other star-forming 
objects ought to be performed, which is now possible thanks to the interactive 
database tool designed by Desert et al. (2006a) and Desert et al. (2006b). 

6.5 Perspectives 

One of the most important complication in the studies of absorption-line spec- 
tra of galaxies is the presence of unseen components, possibly saturated. This 
caveat is usually solved by the use of the weakest observed lines. In that vein, 
the HST/COS instrument enables the sulfur abundance determination in the neu- 
tral phase through the interestingly weak S n A1256 multiplet (Pettini & Lipman 
1995). The oxygen abundance will be measured with a better accuracy than with 
FUSE thanks to the O I A1356 line (Meyer et al. 1998). Furthermore, the exami- 
nation of phosphorus lines will confirm the possible use of phosphorus abundance 
as a metallicity tracer as suggested by Lebouteiller et al. (2005). This is particu- 
larly important for the analysis of the FUSE database since O I lines are generally 
saturated or close to saturation. Finally, it must be stressed that the ideal case 
study for deriving abundances in the neutral gas of BCDs is to observe the ISM 
toward a quasar line of sight (Kunth & Sargent 1986). The method was success- 
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fully used in the dwarf spiral galaxy SBS 1543+593 by Bowen et al. (2005) but 
unfortunately, no alignment BCD-quasar has been observed so far. The study of 
such an alignment would cancel the selection effect in FUSE spectra to probe gas 
toward star-forming regions within the galaxy. 

As far as chemical models are concerned, Recchi et al. (2004) examined the 
case of I Zw 18 for which they found that abundances in the H I medium (defined 
by a temperature lower than 7 000K) are similar to those in the H II gas. Hence, 
although some physics is still missing in the models (radiative transfer, clumpy 
ISM, dust phase), it seems that FUSE results still cannot be reproduced. 

7 Conclusion 

The lower metallicity observed in the H I gas of BCDs as compared to the ionized 
gas of their H n regions remains a puzzle. We have discussed many effects that 
may affect the determination of the neutral gas abundances such as ionization 
corrections, depletion on dust grains, stellar contamination, and line saturation. 
The undcrabundancc of a factor 10 found in the H I region of several BCDs is 
however so dramatic that it seems difficult all these effects can make us miss the 
detection of 90% of the metals. 

It is possible to reproduce the abundance discontinuity in the 2 phases assum- 
ing that only a fraction (typically a few percents) of the metals released by the 
current starburst episode mix within the ionized gas, while most of them mix with 
the neutral envelope. The considerable mass of H I results in only a slight metal- 
licity increase, even with most of the metals mixing within. The main parameter 
therefore concerns the ratio of ionized-to-neutral gas, in the sense that a galaxy 
with a large fraction of its total gas mass ionized is more likely to show similar 
abundances in the 2 phases. Note that his hypothesis is simple speculation at this 
stage and it has to be tested against dispersion spatial scales and against mixing 
timescales. 

It is essential to confirm the existence of a metallicity floor. The lowest metal- 
licities observed in the neutral gas of BCDs agree with the metallicity in I Zw 18 
— the most metal-poor star- forming galaxy known — , suggesting a minimal en- 
richment regardless of the ISM phase. Such a minimal enrichment seems to be in 
agreement with the IGM metallicity around a redshift ~ 2, so that BCDs might 
have formed only "recently" (< lOGyrs ago). 
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